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ABSTRACT 

A  brief  sumraary  of  the  literature  pertaining  to  the  theory  of 
nueues  is  presented  to  acquaint  the  reader  with  the  mathematical 
techniques  avaUahle. 

A  nuantitatire  analysis  of  the  landing  process  around  an  air- 
craft carrier  is  made  vdth  a  view  towards  a  more  efficient  utilisa- 
tion of  the  landing  platform.  The  theory  of  queues  supplies  the 
mathematical  techniques,  thus  making  it  possible  to  codb  to  gripe 
analytically  with  the  problem  of  reducing  the  landing  time. 

The  problem  has  been  idealised  severely  but  is  applicable  under 
certain  conditions  to  the  actual  landing  process.  The  results  may 
be  used  for  investigating  the  distribution  of  the  laeans  of  (a)  the 
aircraft  waiting  time,  (b)  the  time  wasted  by  the  aircraft  carrier 
during  the  landing  operations,  and  (c)  the  number  of  aircraft  waiting 
provided  we  know  (or  assume  we  know): 

(i)  the  minijaum  safe  landing  interval, 

(ii)  the  distribution  of  the  landing  times  (landing  rate), 

(iii)  the  distribution  of  the  time  intervals  between  aircraft 
entering  the  landing  pattern  (departing  rate). 
If  we  have  not  assiuned  the  above,  this  information  most  be  based  on 
data  obtained  from  observation  of  actual  landing  ofperations. 

A  method  is  also  indicated  whereby  the  problem  may  be  analjrTWd 
by  a  quasi-eraperical  technique.  This  is  a  well  known  substitute  for 
an  analjrtical  treatment  if  the  matheiiBtics  become  intractable. 
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niEFACE 

With  the  adTBiit  of  the  jet  age,  terminal  congestion  and  re  stilting 
aircraft  delays  are  of  Increasing  in5)ortance  in  all  branches  of  avia- 
tion. The  ecientifio  approach  provides  a  laeans  for  predicting  how 
these  delays  vflll  be  affected  by  varying  conditions,  and  so  makes  possi- 
ble an  economic  valuation  of  efforts  directed  towards  reducing  these 
delays  • 

Operational  research  methods  have  a  threefold  application  in  this 
connection:   (a)  in  ascertaining  the  wide  field  of  practical  problems 
susceptible  to  treatment;  (b)  in  obtaining  the  basic  data  needed  to 
enable  the  delay  problems  to  be  specified;  (c)  in  selecting  approp- 
riate solutions  i 

oueueing  theory  supplies  certain  mathematical  techniqvies  Y»hereby 
the  operaticais  analyst  can  come  to  grips  analytically  with  the  problem 
of  aircraft  delays.  Even  where  mathematics  become  intractable,  there 
are  good  prospects  of  empirical  solutions  by  means  of  electronic 
calculators . 

The  purpose  of  this  paper  is  to  apply  queueing  theory  to  a  prob- 
lem of  sane  in?)ortance  in  naval  aviation,  carrier  aircraft  landing  de- 
lays. It  follows  that  the  results  may  allow  a  more  efficient  utHi- 
eation  of  the  landing  platform.  Two  things  are  of  particular  interest 
in  the  general  queueing  problem,  the  sise  of  the  queue  (stack)  and  the 
waiting  tiras  of  the  customer  (aircraft  delay). 

Qiapter  1  gives  an  analysis  of  the  general  queueing  problem  and 
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certain  technical  implications  v*iich,   while  they  are  not  strictly 
intef^al  to  the  problem,   throw  considerable  light  on  the  report.     Chap- 
ter 11  presents  a  detailed  analysis  of  one  aspect  of  the  carrier  air- 
craft problem,   a  svmmiary  of  the  reffults,   tlie  conclusions,    and  further 
considerations  of  statistical  implications. 

This  study  was  undertaken  at  the  T^iited  States  Naval  lostgraduate 
School  during  the  latter  hrJLf  of  the  academic  year  1953-54.     "-he  writer 
is  deeply  indebted  to  Frofessor  Charles  C.  Torrance  whose  constant 
criticism  and  encouragement  were  invaluable.     To  him  go  my  foreraoat 
thanks.     The  ivriter  is  also  indebted  to  Professor  A.  Boyd  Ifewbom  for 
reading  the  paper  and  for  making  valuable  suggestions  and  remarks. 
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CH.\PT^  I 
OUEUEING  THEORY 

IJL     Introdxiction  to  the  OxMueing  Frohlem. 

The  Bt\jdy  of  ouexies  has  been  of  interest  to  mathematicians  and 
engineers  for  the  past  forty  years.     <^ueueing  occurs  viien  demands  for 
service  are  greater  than  availability  of  service,   and  it  becomes  neces- 
sary to  postpone  the  demands  by  a  system  of  queueing  or  marshalling. 
'Vhlle  queueing  is  a  natural  process,  marshalling  is  an  attenqrt.  to  do 
something  about  it.     ^ueueing,   however,    is  essentially  a  tenqxjrary 
jiienomenon,   as  otherwise  the  queue  would  grow  indefinitely  in  sise. 
Thus  one  of  the  first  problems  in  any  queueing  situation  is  to  determine 
iirtiether  or  not  the  demand  wHl  outstrip  the  service  mechanism,     "/hen 
this  does  not  happen,   so  that  the  queue  is  continually  returning  to 
aero  size,   then  the  situation  is  said  to  be  in  equilibrium,   or,   in  statis- 
tical terminology,   the  process  is  stationary.     Vftien  the  demand  equals 
the  service  capabilities  the  situation  becomes  very  delicately  balanced; 
a  slight  reduction  in  demand  yields  a  sensible  queueing  system  in  equi- 
librium, vAilLe  a  slight  increase  results  in  an  ever  lengthening  oueue. 

The  theory  of  queues  has  a  special  appeal  for  the  scientist  interes- 
ted in  stochastic  processes,   for  \inder  certain  conditions  it  provides 
an  exanqJLe  vidiich  is  both  stationary  and  non-Miarkovian .     Thus  one  may 
eocpect  in  studying  the  oueue  to  gain  Insight  into  other  stochastic 
phenomena,   and  to  acnulre  a  valuable  facility  in  the  handling  of  the 
relevant  techniques.     In  addition,  the  theory  has  an  astonishing  range 
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of  applicatlona.  A  crirrent  limitation  of  the  theory  is  that  the  num- 
ber of  papers  that  deal  vdth  the  theoretical  side  of  the  subject  is 
small.  And  unfortunately  those  that  exist  are  rather  hea-ry  reading, 
1*2  '"uausing  system. 

The  siii9iLe  nueuelng  system  is  concerned  not  so  BBxch  with  the 
matching  of  supply  and  demand,  but  with  a  sorlos  of  customers  demand- 
ing service  at  a  single  oo\mter  and  waiting  in  turn  to  be  serred.  Cus- 
tomer is  here  used  in  a  technical  sense;  in  the  present  paper  it  should 
be  eruated  with  aircraft.  To  complete  the  specification  (Kendall  2  )  one 
must  gire  a  careful  aoeount  of: 

(a)  the  input-process, 

(b)  the  queue-discipline, 

(c)  the  servlce-meohanism, 
1,2,1  Input-process, 

The  input  is  random  whenever  the  chance  or  probability  that  a 
customer  will  arrive  at  any  instant  of  time  Is  the  same  for  all  instants 
past,  present,  and  fixture,  and  also  is  unaffected  by  any  past  or  future 
arrivals.  The  random  input  is  often  called  a  "polsson  input"  because 
of  the  use  of  the  Foisson  distribution  to  describe  the  arrivals.  A  general 
independent  input  (Llndley  3  )  is  one  where  the  intervals  between  the 
arrivals  of  successive  customers  are  independent  with  identical  pro- 
bability distributions,  biit  this  distribution  is  arbitrary.  The  input 
is  regular  when  the  customers  arrive  at  a  constant  rate.  The  input  is 
correlated  whenever  the  Intervals  between  arrivals  are  not  independent, 
1,2,2  ^ueue-disclpline. 

The  rueue-disclpline  is  the  rule  or  moral  code  determining  the 
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manner  in  vrtilch  the  customers  form  up  into  a  qutua,   and  the  manner  in 
■which  they  behave  v/hile  smiting.     In  the  simplest  case  they  line  up  be- 
fore the  single  counter  and  await  their  respective  turns.     The  theory 
of  multiple  queues  or  many  servers  seems,   «Dccept  under  simplifying  as- 
sun^tions  v;hich  do  not  always  correspond  to  reality,   to  be  a  problem  of 
considerable  difficulty  and  vdJl  not  be  considered  here, 
1.2*3  Service-mechanism. 

The  service^aechanism  may  be  described  as  the  mode  of  service.  In 
this  connection  the  service  time  la  of  basic  Importance  and  is  the  tii» 
that  elapses  while  a  particular  customer  is  being  served. 

1.3  Traffic  density. 

A  parameter  of  special  importance  is  the  traffic  intensity  or  den- 
sity €  t   defined  as  the  ratio  of  the  average  arrival  rate  to  the  average 
service  rate.  This  may  also  be  expressed  as  the  ratio  of  the  mean  ser- 
vice time  to  the  mean  arrival  interval, 

1.4  I/arkovian  properties. 

Before  the  ircthecHtical  theory  is  presented,  an  interesting  quali- 
tative feature  of  the  problem  should  be  commented  on. 

Suppose  that  the  state  of  a  stochastic  system  at  time  t  is  described 
by  a  random  function  X(t)j  the  stochastic  process  is  then  said  to  be  of 
Markov  type  (Doob  9)  i-f  a  kno^iledge  of  the  present  value  of  X(t)  makes 
all  information  about  its  past  history  irrelevant  to  a  prediction  of  its 
future  behaviour.  Such  a  process  may  become  non-Iiiarkcrvian  if  part  of 
the  information  contained  in  X(t)  is  suppressed,  for  a  kno;tLedge  of  the 
previous  behaviour  may  enable  some  of  the  suppressed  information  about 
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the  preaent  state  to  be  recorered.  The  oueuelng  process  vrith  a  Foia- 
sonian  input  and  a  single  server  is  Iferkovlan  if  the  present  state  of 
the  process  is  described  by  the  pair  of  random  variables  q  and  v  ,  >rhere 
q  is  the  instantaneous  qu«ue  size  and  v  is  the  expended  service  time 
of  the  custouisr  at  the  head  of  the  que\»;  in  general  it  ceases  to  be 
llarkovian  if  the  state  of  the  process  is  measured  by  the  queue  size 
alone  (Kendall  2  ).  The  only  exception  to  this  statement  occxrrs  when 
the  seirvice  time  has  a  negative  exponential  distribution, 
1.5  Poisson  input, 

Kendall  (  2  )  gives  a  general  review  of  scnue  points  In  congestion 
theory,  A  presentation  is  made  of  the  Fdlaczek  "equlLibri\m"  theory 
for  the  single  counter  queue  fed  by  an  input  of  the  Foisson  type  and 
associated  ivith  a  general  service  tims  distributicai.  It  is  pointed  out 
that  although  the  stochastic  processes  describing  the  fluctuations  in 
queue  sise  is  not  in  general  Markovian,  it  is  possible  to  work  instead 
with  an  enumerable  Iferkov  chain  if  attention  is  directed  to  the  epochs 
at  which  individual  customers  depart.  The  ergodic  properties  of  the 
chain  are  investigated  with  the  aid  of  Feller  »8  (  S  )  theory  of  re- 
cxirrent  events. 

The  simplest  hypothesis  about  the  input  is  one  i-rtiich  states  that 
the  customers  arrive  "at  randcaa",  the  number  of  arrivals  in  tijne  T 
being  a  Foisson  variable  of  expectation  T/a,  where  a  is  the  mean 
arrival  interval*  The  time  interval  t  between  two  consecutive  arrivals 
will  then  have  the  negative  exponential  distribution 

6F  it)  '-  ^         dl.    ,     o  <  t  <  oo  .  i-i 

This  is  seen  if  one  recalls  that  dF(t)  is  simply  the  probability  of 
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no  arrivals  in  a  tine  interval  t.  The  successive  t-variables  will  be 
statistically  independent.  This  is  the  hypothesis  adopted  by  Kendall, 
1«5J.  lisan  waiting  time  ^vith  general  service  time  distribution. 

Let  the  queue  be  of  the  type  with  a  single  server  fed  by  a  Tois- 
sonian  input  and  let  the  traffic  intensity  €  be  less  than  unity,  so 
that  the  system  is  not  saturated.  Let  q  be  the  size  of  the  ^ueue  which 
the  departing  cxiatoiner  leaves  behind  him,  not  including  himself,  but 
including  the  next  person  to  be  served  (n  may  be  aero),  "^he  next  per- 
son to  be  served  has  a  service  time  v  r']](v)  -  bj  ,  say,  and  suppose 
that  during  this  time  r  new  customers  arrive.  Then  conditionally  r  is 
a  loisson  variable  of  meem  value  v/al  and  v  has  the  service  time  dis- 
tribution. Let  '^'  be  the  siae  of  the  ^ueue  which  the  next  person  leaves 
behind  him. 

i^sun«  that  statistical  equilibrium  exists.  Then  the  random  var- 
iables q  and  q»  have  the  sane  marginal  distribution,  and  in  particular 
their  msan  and  mean  square  values  must  be  equal  and  finite: 

E(%)  ••  rC%'J  <  '>^  3  E  (%M  »  E  (c^'")  <  oo  .      1.2 

The  variables  q  and  q»  are  related  by  the  fundamental  fornula 
^'  .  mux   (V';Q)  -+-  ^  ;  1-3 

or  in  terms  of  the  i  notation, 

%'--    ^-'    ^^    ^^^  1.4 

where  J  -  S  (  %)  is  ''•ro  for  all  non-zero  q,  and  ^  (O)  :  1.  It  is 
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important  to  note  the  foUovdng  conaenuenceo  of  the  definition  of  the 
function   ^  (  %): 

^%  /       and    {c^)U'i  )    --  V-  1-5 

On  forming  the  expectations  of  both  sides  of  Eq,  1.4,  it  will  be 
found  that 


Ei^)-    \'E(n.)    .    \'^A       .   /-  6  . 
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This  is  the  chance  that  i   is  non-zero;  that  is,  it  is  the  probability 
that  an  incoming  customer  ;dJLl  not  have  to  wait.  On  squaring  both  sides 
of  Eq,  IJ^,  and  making  use  of  Eq,  1.5,  one  finds  that 

%'*-  %  -  ^V  t'-'^l^-t^t^^-:^!  (rA.-t)  .  1.7 

Take  expectations  and  note  that  r  is  independent  of  q  and  ^  .  It 
follows  that  ^ 


2  \  \-EiK)\ 


1.8 


Now  suppose  that  a  departing  customer  leaves  q  custcaners  behind 
him,  and  let  his  own  waiting  time  and  service  time  be  respectively  w  and 
V,  Then  q  is  the  nxunber  of  arrivals  in  a  total  tljoB  w  ♦  v.  Thus, 


E(V)  .     ^(^)^^^^)    , 


1.9 


Using  Eqs,  1*8  and  1.9  one  finds  that 
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This  ia  essentially  Kendall's  restilt.  It  is  convenient  to  express 
the  restLLt  in  the  form  of  the  two  ratios,  for  in  applying  the  results  to 
aircraft  delays,  one  usually  normalizes  time  in  units  of  the  mean  arrival 
interval  or  mean  seirvice  time  (minimum  safe  landing  interval), 

"/hen  the  mean  arrival  interval  and  the  mean  service  time  are  equal 
to  each  other  (€  ■  1),  the  expected  waiting  time  is  infinite,  thus  show- 
ing that  the  intuitive  solution  €  ■  1  is  not  the  best  result. 

From  Eos.  1.10  or  IJLl  one  can  see  that  the  minimum  expected  wait- 
ing time  will  be  obtained  if  the  service  time  distribution  is  of  the 
form  /-  -^v.  '  -Sf-v-fc) 

FCv) -- o  C^<t);  FCv^-i  C^^t)^        1.12 

that  is,  the  service  time  is  constant.  On  the  other  hand,  if  the  service 
time  is  of  the  form 

d  F(v)  =  e    iiL  .   o  <  V  <  «o  lj_3 

b 

the  expected  waiting  time  is  twice  the  minimum  value  (voa^^  »  1). 

It  is  tempting  to  think  that  the  expressions  for  E(q)  and  E(w)  may 

be  valid  for  more  general  input  processes;  but  this  seems  to  be  a  fallacy. 

The  difficulty  is  that  r  and  q  are  no  longer  statistically  independent. 
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Kendall  shows  that  the  ^irkov  chain  aaaociated  with  the  qxaeueijig 
process  is  irreducible  and  the  states  are  aperiodic.  A  further  classifi- 
cation of  the  states  depends  on  the  ralue  of  the  traffic  intensity,  €  • 
\'lt»n   €  is  less  than  unity,  all  states  are  ergodic;  when  €  i«  greater 
than  unity,  all  states  are  transient;  in  the  critical  case  when  €  is 
equal  to  unity,  the  states  are  all  recurrent  and  nuH.  (Feller  9  ) 
1,5,2  AJLrcraft  delays. 

Bowen  and  Fearcey  (  6  )  give  an  elemsntary  analysis  of  the  effect 
of  control  procedures  on  the  flow  of  air  traffic,  while  Fearcey  (  7  ) 
presents  a  more  rigorous  treatnent  of  the  problem  and  derives  the  dis- 
tribution of  delays  under  different  traffic  conditions.  The  hypothesis 
of  random  arrivals  is  used,  so  that  one  is  led  to  the  negative  exponen- 
tial distribution  for  the  arrival  intervals.  See  3q.  IJ.. 

The  assumption  is  made  that  the  minimum  safe  landing  interval 
(service  time)  is  consrtant.  Thus  we  have  Kendall »s  solution  with  a  con- 
stant service  time  in  a  different  notation. 

Fearcey 's  treatment  is  based  on  a  rather  con^xLicated  procedure  of 
successive  integrations  but  once  the  assumption  is  made  that  the  system 
is  in  a  state  of  statistical  equilibrium  the  analysis  most  lead  to  the 
sane  results.  One  feature  of  Interest  is  the  derivation  of  the  proba- 
bility of  successive  delays.  By  successive  integrations  Fearcey  finds 
that  in  general  the  probability  of  n-1  aircraft  being  successively 
delayed  is  (stack  of  n-1  aircraft) 


P(y^)  -.  AC^)T     'e"" 
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""  n-    ,  ,  n. 


""^Ti —  ^ 

Dell  (    5  )  gives  a  somsAvhat  more  general  treatment  by  using  a  moment 
generating;  procedure. 
1,6    General  independent  input. 

The  center  of  interest  of  Lindley's   (   3  )  paper  ±a  the  waiting  times 
of  the  customers.     A  relation  betvreen  the  waiting  times  of  successive 
customers  is  found,   and  this  makes  it  possible  to  calculate  the  wait- 
ing time  distribution  of  any  customer. 

This  relation  is  also  used  to  obtain  necessary  and  sufficient  con- 
ditions for  the  existence  of  a  stationary  (e  uilibrium)  state.  The  paper 
concludes  ^vith  an  investigation  into  the  artationary  state  in  some  spec- 
ial cases;  particularly  in  the  case  when  the  cvistomers  anrive  at  regular 
intervals. 
1.6  J.  ,/aiting  time  distribution. 

The  -^lUeue  here  is  again  of  the  type  with  a  single  server  but  with 
general  service  time  and  input  distributions;  the  intervals  between 
successive  arrivals  may  have  any  suitable  distribution. 

Let  tj.  be  the  time  interval  betv/een  the  arrival  of  the  r       and 
p  •♦■  1       customer,  and  let  s     be  the  service  time  of  the  r       customer. 

The  following  assumptions  are  made: 
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I4  The  t^  are  independent  random  variables  vrith   Identical 

probability  distributions  and  the  Man,  E(t  ),  is  finite: 

r 

oO 

^^^r.)      =     J    t^      pCt^]      dt^      <      CO 

o 

2,     The  8     are  independent  random  variables  vdth  identical 
r 

probability  diatribtttions  and  the  mean,   E(8  ),   ie  finite. 

Thus  the  two  sets  of  random  variables  [s    ]    and    [t  ]       (r  -  1,2, .,.) 
are  atatistically  independent. 

Let  w    be  the  waiting  time  of  the  r       customer.     It  is  evident 
that  the  following  relations  exist  (Fig.  1«1): 


n.-».i 


n4l 


5  O 


n.4'1 


Fig.  1,1     V/aiting  Tijns  Relations 
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(t,  <  »r  *  'r  ) 
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r  ♦  1 


(tj.  ^  1^  ♦  8y;  ■•e  ty) 
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Let  Up  ■  8  -  t  -  The  u,  are  independent  random  variables 
(by  virtue  of  the  proceding  assumptions)  and  have  identical  probability 
distributions;  E(  |  Uy  |  )  is  finite  and  u^  is  independent  of  w^..  Then 


r  +  1 


-   Wj.  ♦Uy 


•  0 


(Wy  ♦  Uj.  >  0), 
(Wy  ♦u^  10). 
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Now  if  ym  aussxziaa  that  th«  first  cuatomer  does  not  hscre  to  wait 
(w,    -  0)  we  have  the  following  relations: 


and  the  unexpected  fact  develops  that  the  ?(aitljDg  t lines  depend  on  only 
the  difference  between  the  service  time  and  interval  time.     If  we  define 
F  (x)  as  the  probability  that  w^  i  x,   and  G^(x)  as  the  probability  that 

T  T  T 

u    ^  X,  Tf»  can  sll^tly  modify  Lindley»s  procedure  and  find  the  dlatri- 
r 

art 
butlon  of  Tfraiting  times  of  the  r  ♦  1       customer  by  using  Eq.  I.IS.     The 

u    are  independent  statistically.     Thus  the  distribution  function  of 

Fj,  ♦  1  is  the  convolution  of  u    ♦  u^^  •»■ ♦  il    (Doob  9)     and 

pK„  ^  ^)  ■■    i.t  G.  Cu.)  .  .  .  .  XJ^  (  k-  u, u^.,)  dG^.,(V,\^5 

Using  the  definition  of  Gj.(x),  we  have 

Eq.  1.19  establiahes  a  recursive  relation  such  that  the  individ- 
ual waitinn;  time  distributions  may  be  found.  To  find  the  final  behavior 

of  the  cueue,  we  must  find  the  limit  of  F  ^,  (x)  as  the  number  of  cus- 

r*l 

tomers  continually  increase.     Let  E_  be  the  event: 
*  r 

Now  the  sequence  of  events  [lilpl  decreases  to  a  limit  E  as  r  approaches 
infinity,  that  is: 
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by  a  theorem  of  probability  maasiire.     In  the  limit,   then, 


1.24 


This  integral  equation  is  difficult  to  solve,   but  it  is  possible  to  do 
so  once  G(u)  is  specified  and  certain  conditions  are  satisfied     (Paley  and 
Weintr  IL).     The  solution  of  either  integral  equation  gives  the  dis- 
tribution function  of  the  waiting  time,  that  is,  the  probability  of 
waiting  a  time  not  aore  than  x.     The  distribution  function  for  any  par- 
ticular customer  can  be  calculated  by  performing  the  necessary  inte- 
grations; but  this  is  also  difficult  to  do. 

Lindley  uses  Feller's  theory  of  recurrent  events  and  the  strong 
law  of  large  numbers  to  investigate  the  behavior  of  F(x), 

The  event  of  "not  having  to  wait",   called  event  A,   is  a  recurrent 
_  event,   since  the  queueing  process  may  be  considered  to  make  a  fresh  start 
whenever  a  customer  does  not  have  to  wait,  and  its  previous  history  is 
quite  irrelevant  to  subsequent  developments* 

To  investigate  the  behavior  of  F(x)  one  most  examine  the  relation 

E(u).E(s)   _   ECO  ,.,5 

as  fbllcnra: 

E(u)  >  0:     The  strong  law  shows  that  the  event  A  is  transient  and 
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that  with  probability  1  there  vdll  be  some  point  in  the  queuelng  process 
irtiere  every  subsequent  cuatomer  will  have  to  wait, 

E(u)  <  0:  The  strong  law  shows  that  event  A  is  certain  with  a  finite 
mean  recurrence  tiioe.     It  should  be  clear  that  it  is  not  periodic. 

E(u)  -  0:     The  event  A  is  certain,  but  the  mean  recurrence  time  is 
infinite.     Hence  it  is  certain  that  there  will  always  be  another  occasion 
on  which  a  customer  does  not  have  to  wait  but  the  probability  of  any 
particxilar  customer  not  having  to  wait  tends  to  zero. 

Thus  a  necessary  and  sufficient  condition  that  the  waiting  time  dis- 
tribution function  tends  to  a  non-degenerate  limit,  as  the  number  of  cus- 
toaaers  increases,   is  that  E(u)  <  0  or  u  -  0  certainly, 

Lindley  solves  the  integral  equation  when  the  arrivals  are  random 
that  is, 

•    dGCt)   =>>   «  dt     ;    ^'^^^--t    >  1.26 

and  also  for 

d  &(^)   =    x\  e  '^^d^  •     El^)  =  X  .  1,27 


These  distributions  are  the  Type  HI  curves  of  Pearscxi  ccaononly  called 
the  garma  di«tributicns  (Mood    10),  and  are  of  the  form 

ri  !       "  ^  '^  X 

with  the  given  parameters  n  and  X  • 
And  as  expected,  the  probability  of  not  having  to  wait  is 

1-  111}     •.\'  e  . 

ECt) 
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The  •quation  to  be  solved  for  regular  arrivals,  where  the  arrival 
interval  is  the  unit  on  the  time  scaler  is 

,  F'^)  =  i^f,  '^^^-^^=*^.  (^^0-  1,30 

vdth  the  service  time  distribution  as  given  by  Eq,  1,2S, 

The  solution  is  difficxiLt  to  develop,  but  many  may  be  expressed  in  the 
form 


1,32 


n  +  i               3.x 

F  U)  =   1-   Z:   Ci  e  ^ 

) 

C  .^  :    C  .        and 

1                    '^■^  *      ^  . 

o 

^■"'       '■■'    (x^^)'^^' 

K^ 

^  ,  =  Z-,     9  n  d 

^      =    e 

The  usual  form  of  service  time  distribution  is  one  whose  frequency 
function  is  zero  at  s  ■  0«  increases  to  a  maximum  and  then  decreases, 
with  a  long  tail,  to  zero  as  s  approaches  infinity.  In  seme  cases  it 
may  decrease  steadily  f^m  some  non-zero  value  at  s  -  0«  To  a  good 
approximation  these  functions  are  represented  by  the  gamma  distribu- 
tions (n  -  0,  1), 

The  results  of  the  calculations  show  that  regular  arrivals  are 
roughly  twice  as  efficient  as  random  arrivals,  i,e«,_thi_j!robability  of 
not  having  to  wait  for  regular  arrivals  is  roughly  twice  as  great  as 
the  random  arrivals.  Lindley's  calculations  for  a  few  values  of  ^  are 
listed  in  Table  1  to  show  the  nature  of  the  results  (units  ai^  arbitrary). 


Probability  of  not 

having  to  wait 

n  -  0 

Regular 
Arrivals 

6 

Random 
/iTrivals 

0.830 

0.314 

oa67 

0.715 

0.511 

0.286 

0.635 

0.642 

0.375 

0.500 

0.797 

0.500 

n  -1 

^ 

0.830 

0.386 

0.167 

0.715 

0.608 

0.286 

0.635 

0.743 

0.375 

0.500 

0.884 

0.500 

Ifean  waiting  time 

Regular      Random 
Arrivals     .^jrrivals 


1.823 

4a66 

0.684 

1.786 

0.349 

1.042 

0.128 

0.500 

0.834 

3.124 

0.285 

1.339 

0^.32 

0.781 

0.039 

0.375 

Table  1 


1.7     Correlated  input. 

In  view  of  the  difficulties  encountered  vdth  the  random  and  inde- 
pendent inputs,   it  should  be  readily  apparent  that  the  inclusion  of  con- 
ditional probabilities  in  the  input  process  vdll  make  the  protlem  even 
more  difficult,      ^ji  mi^t  be  expected  then,  this  problem  has  not  been 
solved  analytically. 

The  use  of  the  Poisson  input  in  air  traffic  leads  to  a  complete 
analytical  solution  of  the  stack-delay  problem,   but  unfortunately  once 
aircraft  are  in  the  air,  the  input  is  not  even  approximately  random, 
but  has  the  form  of  a  regular  pattern  vdth  large  errors.     It  is  virtually 
in^ossible  to  find  the  actual  distribution  vdth  any  accuracy,   as  the 
amount  of  data  required  for  an  unbiassed  sample  is  prc^bitive.     The  key 
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to  the  difficulty  lias  in  the  fact  that  aircraft  in  the  air  proceed  on  a 
predetermined  schedule.  This  detailed  information  in  the  schedule  should 
be  regarded  as  a  datum  of  the  problem  and  used  in  its  solution* 

The  first  significant  contribution  of  the  electronic  calculator  to 
the  que\jeing  problem  has  been  its  use  in  the  study  of  the  flow  of 
scheduled  air  traffic.  The  Research  Laboratory  of  Electronics  at  the 
Massachusetts  Institute  of  Technology  {  U  )   are  the  pioneers  in  this 
endeavor.  They  undertook  a  program  to  determine,  as  far  as  practicable, 
the  quantitative  relation  betiroen  the  degree  of  control  of  aircraft 
enroute  to  an  airport,  and  the  resulting  congestion  and  delay. 

Actually  there  are  t-wo  possible  problems  involved  in  the  control 
aspects  of  air  traffic «  One  is  the  maintenance  of  a  pre-assigned  schedule, 
and  the  other  is  the  reduction  in  the  congestion  arovmd  the  terminal  or 
landing  strip.  It  should  be  noted  that  this  is  an  unusual  situation  in 
that  the  aims  of  the  customer  and  the  server  are  the  same:   (a)  to  ensure 
safety  in  the  airj  (b)  to  allo«f  the  most  efficient  utilization  of  the 
landing  platformj  (c)  to  reduce  the  total  time  by  avoiding  long  delays 
at  the  terminal.  This  is  not  always  true  in  the  general  queueing  problem. 

In  the  MTT  report,  IBM  punched-card  machines  were  employed  to  make 
a  theoretical  analysis  of  the  resulting  congestion  at  a  single  landing 
strip  vhen  aircraft  are  scheduled  to  arrive  there  in  some  proper  se- 
quence, but  fail  to  meet  such  schedules  according  to  certain  simple 
deviaticm  statistics  (rectangular,  triangxalar,  and  parabolic).  These 
deviations  statistics  are  assumed  to  have  a  finite  spread,  that  is, 
deviations  in  excess  of  a  specific  amount  have  sero  probability.  This 
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is  within  the  bovmds  of  reason  from  the  pilot's  viewpoint.  Numerical 
results  for  the  distribution  of  the  resulting  staok  and  total  delays 
are  presented.  The  former  are  compared  with  those  which  arise  from 
the  more  random  Foisson  arrival  distribxition.  The  amounts  of  conges- 
tion and  stack  delay  are  foxmd  to  be  considerably  smaller  under  the 
new  conditions,  particularly  when  the  traffic  is  heavy.  A  significant 
range  of  the  parameters  has  been  covered  by  the  numerical  analysis » 
supported  by  strictly  analytical  methods  wherever  possible.  Thus  there 
is  establidied  a  quantitative  relation  between  the  time-keeping  errors 
of  the  aircraft  enroute  to  an  airport  and  the  restilting  terminal  con- 
gestion and  delay. 

The  density  of  the  traffic  flow  to  an  airport  is  described  by  the 
usual  traffic  parameter  €  ,  defined  as  the  ratio  of  the  average  arrival 
rate  at  the  airport  (number  of  planes  par  \2nit  time)  to  the  maximum 
allowable  acceptance  rate.  This  is  also  the  ratio  of  the  landing  in- 
terval (service  tirns)  to  the  arrival  interval. 

The  authors  found  that  it  was  convenient  to  use  time  in  discrete 
units  rather  than  on  a  continuous  basis.  The  minimum  safe  landing 
interval,  t^,  is  chosen  as  the  unit  of  time,  and  forms  the  base  for 
both  the  normalization  smd  quantization  of  all  other  times  appeeuring 
in  the  problem. 

A  sample  size  of  1,000  planes  was  selected  for  use  in  the  numerical 
analysis.  In  the  programning  procedure,  each  plane  in  the  sample  is 
represented  by  a  card,  which  is  identified  by  a  number  J,  with  J  greater 
than  or  equal  to  1  and  less  than  or  equal  to  1000.  Thus  J  defines  the 
scheduled  arrival  order.  The  scheduled  arrival  times  t .  are  then  as- 
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signed  using  a  random  nuraber  table. 

Between  its  scheduled  take-off  and  actual  arriyal  times  the  ali^raft 
is  subject  to  a  delay  Ta,   which  is  assigned  (punched  onto  the  j   card) 
according  to  the  enroute  distribution  used.  Thus  plane  J  actually  ar- 
rires  at  time  p^  -  t.  ♦  r.,  which  number  is  also  entered  on  the  J 
card.  It  is  possible  now  for  more  than  one  plane  to  arrire  at  the 
same  time  and  for  them  to  arriTS  in  a  sequence  idiich  differs  frora  that 
scheduled.  This  leads  to  the  formation  of  a  stack  of  aircraft,  and 
consequent  delays  while  waiting  in  the  stack.  It  is  then  assumed  that 
planes  will  be  handled  at  the  terminal  in  order  of  actual  arrival,  with 
the  prorlsion  that  simultaneous  arrivals  will  be  considered  as  ordered 
according  to  their  originally  scheduled  sequence. 

After  the  preliminary  details  were  completed,  i.e.,  the  t  and  r 
tables  were  prejared,  and  the  initial  conditions  were  chosen,  the  IBM 
machines  were  programmed  in  detail  for  the  specific  cases,  i.e.,  for 
each  deriation  distribution,  a  spread  S,  and  traffic  parameter  €  • 

The  principal  results  which  caiB  directly  from  the  machines  were: 

(a)  Frequency  distributions  for  the  stack  delays. 

(b)  Average  stack  delay  for  each  nm. 

(c)  Cumulative  frequency  distributions  for  the  stack  delays. 

(d)  Frs'-uency  distributions  for  the  total  timekeeping  errors 
(stack  delay  plus  enroute  delay). 

(e)  Average  timekeeping  error  for  each  run. 

{ty    Cumulative  fre-usncy  distributions  for  the  total  timekeeping 
error. 
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I'dthout  loss  of  generality,  it  vras  assumed  that  all  variations 
were  delays  (simple  change  of  scale). 

The  results  form  a  considerable  source  of  data  pertinent  to  a 
traffic  of  pj^perly  scheduled  aircraft.  The  data  are  given  in  detail 
in  the  form  of  curves  of  stack-delay  probabilities,  prcglressive-delay 
probabilities,  and  average  delays,  for  the  three  types  of  enroute-de- 
viation  distributions  considered. 

The  curves  give  additional  information  regarding  the  effects  of 
variation  of  the  parameters.  In  particular,  only  relatively  small 
differences  in  the  results  are  caused  by  rather  moderate  variations 
In  the  shape  of  the  enroute-deviation  distributioiis .  The  relatively 
small  changes  in  the  stack-delay  distribution  caused  by  making  moderate 
changes  in  the  shape  of  the  enroute-deviation  for  given  S  made  it 
possible  to  average  together  the  IBM  results  for  all  th3?ee  basic  shapes. 
The  average  curves  consequently  give  stack-delay  probabilities  for  an 
average  shape  enroute«<leviaticn  distribution,  iriiioh  is  sort  of  cross 
between  the  rectangular,  triangular,  and  parabolic  shapes. 

The  effects  of  changes  in  the  traffic  parameter  6  are  practically 
the  same  for  all  three  types  of  distribution.  For  values  of  €  close 
to  unity,  the  average  delay  changes  fairly  rapidly  vdth  €  ,  wrtdle  cor- 
responding changes  in  €  ,  vdien  it  has  smaller  values,  cause  little 
change  in  the  average  delay.  The  stack-delay  distribution  is  also  de> 
pexxlent  in  a  similar  way  upon  the  traffic  parameter.  %en  6  is  close  to 
unity,  the  most  probable  stack  delay  is  almost  equal  to  its  average  value, 
but  as  €  decreases,  the  most  probable  delay  rapidly  approaches  zero 
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with  a  small  average  delay.  The  remaining  rariahle  parameter,  the  apread 
S,   also  has  a  decreasing  effect  with  decreasing  ralues  of  €  « 
l.d  Tabiilation  of  eases  covered. 

The  principal  oases  covered  may  be  tabxiLated  briefly  as  follows: 

(a)  Poisson  input  with  general  service  tijne  distribution  (Kendall  2  ). 
Derives  expressions  for  the  expected  waiting  time  and  the  ex- 
pected qxisiia  siae. 

(b)  General  Independent  Input  i*e,,  general  service  time  and  input 
distributions  ineluding  regular  arrivals  (Llndley  3  )• 
Derives  the  waiting  time  distributi(m  of  any  cusrtoner  and  the 
distribution  function  of  the  waiting  time  in  the  statioiary 
case. 

(c)  Correlated  input  i.e.,  scheduled  air  traffic  (Adler  and  Pricker 
4   ).  Use  electronic  calculator  (IM)  to  find  the  frequency 
distxrlbutions  for  the  stack  delays,  cumulative  frequency  dis- 
tributions for  the  stack  delays,  and  the  average  stack  delay 
for  each  run.  The  information  is  presented  in  the  form  of 
useable  graphs. 

1,9  ImpQ.  icat  ions 

An  interesting  interpretation  of  the  MTT  results  may  be  visualised 
in  the  following  manner.  Consider  the  problem  of  bombing  when  a  large 
number  of  planes  may  be  directed  to  strike  several  targets.  The  usual 
practice,  to  avoid  saturating  the  handling  capacity  of  the  airfield,  is 
to  stagger  the  formations  in  such  a  manner  that  all  of  the  aircraft  will 
not  return  at  the  same  time. 
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For  exarapLe,   consider  a  200  jUaxM  strike  made  up  of  20  groupe  of 
10  jian*s  per  group.     Now  m  hare  no  Idea  of  the  enroute  derlatlon  dis- 
tribution,  but  let  us  use  our  ImaglnaticHi  and  say  that  the  average 
enroute-deriation  distribution  described  in  Sectiosi  1«7  la  applicaKLe. 

Assume  the  expected  fli^t  time  to  target  and  back  to  be  10  hours 
for  each  group.     Let  the  mlnlnwim  safe  landing  Interral  be  approximately 
2  minutes  per  plane  so  that  it  takes  on  the  arerage  20  minutes  to  land 
a  group  (t     ■  20  minutes). 

It  is  desired  that  the  20  groups  arrive  home  in  such  a  manner  that 
the  planned  maximum  stack  is  to  be  2  groups  or  less  (20  planes).     Larger 
stacks  are  not  feasible  due  to  lack  of  control  facilities,  safe  altitude 
separation,   etc. 

For  a  10  hour  fli^t,  errors  of  *    20  raimites  would  not  be  deemed 
too  large;  in  the  notation  of  the  inT  report,  this  vdll  correspond  to  a 
spread  of  i^O  minutes.     If  ire  normallBS  the  spread  in  units  of  t^,   S  -  2. 

If  we  apply  the  results  of  the  JJHH  paper  with  €   ■  1,  that  is,  with 
the  departure  IntervtCL  also  20  minutes  and  with  S  ■  2,  we  find  the  pro- 
bability of  stacks  greater  than  or  equal  to  2  groups  is  approximately 
0.95*     In  other  words,   it  is  almost  certain  that  more  than  2  groupe  will 
overlap  at  the  field.     If  we  reduce  6   to  0.8,  that  is,  make  the  depar- 
ture Interval  25  minutes,  the  probability  that  stacks  greater  than  or 
equal  to  2  groups  will  form  is  reduced  to  0.06.     Thus  it  may  be  seen 
that  queueing  theory,   as  developed  by  electronic  computers  used  in  con- 
Junction  with  good  experimental  data,  has  strlkliig  possibilities. 

This  eacampLe  is  of  course  trivial  in  that  it  is  such  a  restricted 
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problem.  But  enren  ao,  the  results  are  well  worth  further  study.  The 
Polsion  solution  whan  €  is  1  has  been  stressed  enough;  for  example,  the 
stack  approaches  infinity,  and  the  waiting  time  appi»oaches  infinity. 
Thus  there  would  be  no  solution  in  that  event.  IThen  €  is  0,S,  the 
Poisson  solution  gives  for  the  probability  of  2  or  imxre  groupe  in  the 
stack  a  value  of  0,25  apprcodmately.  In  coiaparing   the  Foleaon  solu- 
tion with  the  IHI  solution  we  see  the  two  results  differ  by  a  factor 
of  4>  a  considerable  discrepancy.  The  key  to  the  difficulty  is  the 
use  of  schedule  infonnation  as  a  datum  of  the  problem. 


22 


CHAPTER  n 

C.^IRHILR  .^RCR/vFT  LANDING  DELAYS 
2  J.  Formulation  of  the  problem. 

The  prime  desideratum  in  landing  aircraft  aboard  a  carrier  is 
simply  to  land  the  aircraft  as  quickly  as  possible  vrtdle  observing  the 
necessary  niles  of  safety.  The  time  required  to  land  (and  launch)  the 
aircraft  aboard  the  aircraft  carrier  is  of  basic  importance  due  to  the 
aoctreme  vulnerability  of  the  carrier  during  the  landing  operations.  The 
intuitive  solution  of  the  problem  of  minimizing  the  landing  time  is  to 
minimise  the  landing  intervals  by  feeding  the  aircraft  into  the  landing 
pattern  as  quickly  as  possible  subject  to  the  restraint  of  minimum  safe 
landing  interval.  V/hile  this  is  sometimes  a  satisfactory  solution,  yet 
sometimes  the  delays  involved  become  excessive.  Thus  in  attempting  to 
minimize  the  landing  tirae,  one  must  investigate  the  behavior  of  the  air- 
craft dtiring  the  landing  process.  This  fact  makes  necessary  a  quanti- 
tative analysis  of  the  landing  process,  with  a  view  towards  a  more 
efficient  utilization  of  the  landing  platform.  Thus  it  seems  natural 
to  apply  the  theory  of  queues  in  the  analysis.  Let  us  consider  the 
main  aspects  of  the  situation  under  the  three  headings  introduced  by 
Kendall  to  characterize  a  queue ing  process. 

As  a  preliminary  step,  consider  the  pattern  of  events  around  an 
aircraft  carrier.  All  of  the  aircraft  waiting  to  land  are  placed  in  an 
orbit  above  the  carrier  at  a  safe  altitude  until  the  signal  to  begin 
larriing  is  given.  This  may  be  visxialized  as  a  ready-made  stack  (queue). 
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Vflien  the  signal  to  land  is  giren,  the  aircraft  descend  to  the  laiKling 
pattern  to  oanmanoe  their  approach.  There  is  a  constraint  in  that  onLj 
a  maxinnim  number  of  aircraft  are  allowed  in  the  landing  pattern.  Srery 
time  an  aircraft  lands  it  is  replaced  in  the  landing  pattern  by  one  of 
the  waiting  aircraft.  Thus  one  can  vistialise  the  single  aerver  queue- 
ing  process. 

First,  there  is  the  input  process.  In  our  context  this  is  essen- 
tially the  process  by  vdiich  the  aircraft  depart  the  orbiting  stack  (queue) 
and  feed  into  either  the  landing  pattern  directly  or  Into  another  stack 
(queue)  of  aircraft  waiting  to  enter  the  landing  process.  The  landing 
pattern  may  also  be  visualised  as  a  queue  with  its  associated  pr^- 
blem  of  delays  due  to  wave  offs,  etc.  This  is  a  problem  for  further 
study  and  is  not  considered  here.  Thus  we  may  yisuallBe  two  queues 
with  a  single  server.  It  should  be  noted  that  the  orbiting  stack  (high) 
is  called  a  queue  In  a  gtrictly  limited  sense.  This  artifice  makes  the 
analysis  more  understandable.  The  low  stack  (queue)  is  the  feature  of 
interest  in  that  it  is  caused  by  sod»  perturbation  in  the  landing  process. 
VJ9   shall  assume  that  all  the  aircraft  leare  the  high  stack  at  a  con- 
stant rate  (continuous  flow),  i.e.,  at  a  succession  of  equally  spaoed 
times,  the  difference  between  one  time  and  the  next  being  called  the 
departure  interval.  Thus  the  departure  distribution  is  of  the  form 

where  r  is  the  departure  rate  (constant).  Second,  there  is  the  queue 
discipline.  Mb   shall  require  that  the  aircraft  land  in  the  order  in 

2k 


which  they  arriva  at  the  low  stack  if  one  exists.  Third,  there  ia  the 
service  mechanism,  tihich  is  given  by  the  frequency  distribution  of  the 
landing  times  (serving  time).  In  a  limited  sense  m  shall  consider  the 
landing  pattern  as  a  part  of  the  service  mechanism.  LaiKiing  tirae  is 
here  taken  as  the  time  from  the  cut  signal  given  by  the  landing  signal 
officer  to  the  time  the  aircraft  is  safely  forward  of  the  barriex^,  ''.'}• 
vdll  new  make  the  drastic  simpUf^ring  assumption  that  the  landing  time 
is  constant,  i.e.,  the  landing  distribution  is 

FU)  =  o    (Us)  •    P(s)  .-  I    (j^i  s)^ 

vAiere  1  is  the  landing  rate  (constant).     This  assuioptlon  is  made  for 
two  reasons:     (a)  observational  data  are  lacking,  and  (b)  it  is  desired 
to  make  the  analysis  traetable.     It  seems  rather  natural  to  assume  that 
the  landing  time  distribution  is  one  nhose  fl^quency  function  is  cut  off 
at  the  minijiium  safe  landing  Interval,   increases  to  a  msiximum  near  some 
average  value  of  landing  time,  and  then  decreases  to  another  arbitrary 
cut  off  point  as  the  landing  time  Increases.     Some  consideration  wUl 
be  given  later  to  certain  statiatical  ImjiLicaticoa  pertaining  to  the 
gathering  of  observational  data.     A  second  drastic  assumption  to  be  made 
in  connect icvi  with  the  landing  time  is  that  under  normal  condltlona 
(vrhen  the  normal  laming  rate  is  in  effect)  the  aircraft  are  landed 
aboard  without  any  delay  for  a  certain  length  of  time,   and  then  the 
normal  landing  rate  Is  changed  abn^ly,   and  c(»itinues  at  a  new  reduced 
constant  rate  for  a  period  of  tiias  T^.     The  reason  for  this  reduction 
In  rate  may  be  poor  handling  technique  aboard  the  eairier,   ineocperienced 
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pilots,  etc,     Dxiring  this  period  of  tine  the  low  stack  begins  to  form. 

It  is  noir  easily  seen  iriiy  we  considsr  two  queues  in  series  with  a  single 

server. 

An  in^rtant  parameter  of  the  problem  is  the  number  N  of  planes 

Initially  in  the  high  stack.     Uhder  nonnal  conditions  (ivhen  the  low  stack 

is  empty)  the  time  for  landing  N  aircraft  is 

Ty     .    N/r. 

Thus  the  problem  is  to  inrestigate  the  delays  under  the  two  conditions » 

T   >  T     and  T  t    T^. 
o    r        o    r  • 

Before  proceeding  ^rith  the  analysis  it  should  be  noted  that  the 
fundamental  difference  betTreen  the  general  queueing  problem  and  the  pro- 
blem at  hand  is  that  vre  are  not  Interested  in  the  limiting  case,  that 
is,  in  an  infinitely  long  queueing  process.  Ths  stack  (queue)  assoc- 
iated with  the  carrier  problem  has  a  finite  bound  N  as  a  datum  of  the 
problem,  and  this  must  be  included  in  the  analysis.  Howerer,  as  the 
departure  Interval  approaches  the  landing  time,  we  might  expect  from 
the  work  of  Liudley  that  the  situation  is  somsv^at  oritieal  for  short 
queues  as  viell  as  for  long  ones.  This  is  in  fact  the  case. 

The  notation  to  be  used  is  as  follows: 

X  (t):  landing  rate  aboard  aircraft  carrier  (planes  per  unit  of 
time). 

r  (t):  departure  rate  from  hi^  stack  (planes  per  unit  of  time), 

£  ;  ~  :  normal  traffic  density  parameter  also  equal  to  the  ratio 
of  the  landing  time  to  the  departure  interval, 

a    :  landing  rate  reduction  parameter. 

T    :  total  time  a  is  diminishing  the  landing  rate, 
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T   :  total  time  during  which  nc»i-zero  low  «tack  exiats. 

Ty  :  total  time  for  carrier  to  complete  l«mdings  under  normal 

condltiona . 

N   :  number  of  planes  initially  in  hi^  stack 

S(t):  nuEber  of  planes  in  lew  stack  (low  stack  sisae). 

T   :  total  additional  time  reaulred  for  carrier  to  conmLete 
ao 

landings  due  to  formation  of  low  stack, 
a*   :  aircraft  delay. 
D   :  total  aggregate  delay  (say  in  plane-mimttes)  equal  to  product 

of  average  delay  ( 'T  )  and  the  total  number  of  planes  delajred} 

also  eoual  to  the  product  of  T  and  S(t). 

0 

Note  that  the  traffic  density  €  ia  not  a  function  of  time. 

The  numerical  quantities  of  primary  interest  will  be:   (1)  the  ad- 
ditional time  reqtdred  for  the  carrier  to  con^xLete  landings  due  to  tho 
reduction  of  the  landing  rate;  (2)  the  average  stack  delay  over  the  per- 
iod during  \vhich  a  non-eerc  low  stack  eodsts;  and  (3)  the  maxijmim  stack 
delay  sitffered  by  any  plane. 

The  calculatiana  fall  into  tvvo  major  groups,  defined  by  T  >  T  and 
T  <  Tp.  The  first  case  to  be  considered  vdll  be  T^  ^  T  .  Figs.  2 J., 
2,2,  2,5,  and  2.6  show  the  corresponding  departure  rate,  landing  rate, 
and  stack  height  for  prescribed  values  of  a  and  €  in  the  ranges 

0  <  a  <  £  <1. 

\le  can  assume  without  any  loss  of  generality  that  T  eomnences 
when  t  -  0.  Prior  to  this  time,  the  landijig  process  is  normal,  say 
(  €  t  H.  •  ^^   I  )  and  there  is  no  low  stack.  It  is  not  essential  to 
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the  analysis  whether  the  landing  process  is  in  progress  at  t  -  0,  or 
Just  commencing. 

Now  at  the  beginning  of  time  Tq  the  landing  rate  is  reduced  to  aj^ 
and,  provided  a  <€  ,  a  low  stack  begins  to  form  at  a  rate  (r  -  a  Ji)«  . 

r-a  i(  .  J((€  -a). 
It  is  essential  to  keep  in  mind  that  no  stack  forms  at  all  unless  a  <  e , 
and  this  condition  is  therefore  implicit  throoghoub  the  derivations.  At 

time  T  ,  the  departure  rate  from  the  hi^  stack  beo(»ne8  zero,  and|  in 

r 

ocffisequence  the  stack  no  longer  increases  bat  begins  to  decrease  at  a 

rate  ax  .     .\fter  a  ftorther  period  of  time  (T  -T  ),  the  aircraft  carrier 

0    r 

resumes  its  normal  landing  rate.     As  the  hi^  stack  has  cleared,  the  low 
stack  decreases  vvlth  the  normal  landing  rate  J[  • 
2.2     Stack  siee  S(t)  aixl  stack  time  T.. 
At  time  T,.  the  low  stack  hei^t  is 

S(T«)  i*  eoual  to  S(T)         .     See  Fig.  2J.«     Now  the  stack  begins  to  de- 

'  max* 

crease  at  a  rate  a  jl  ,  and  at  time  T^, 

:  X,  J   ( t  -  c^  '^yr^  )  . 

Now  S(T   )  is  non-negative  if  and  only  if 
o 

L  >  li     • 

a      "*     -r         ) 

'    n. 


2.2 


if  this  is  not  the  case  and 


-5"    ^  rr-       ) 


the  stack  clears  entirely  prior  to  the  end  of  time  T^«  See  Fig.  2«2. 
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Th\2B  yn  have  two  cases  to  Investigate: 


Case  A:     —   ^   H    >   I  .  Caae  B:    Xsl   ^  A-   >  I 


Case  A:  ^  >   TV    >   I 

The  time  taken  to  clear  the  stack  left  at  time  T     Is 

o 


Norr  the  tins  that  a  nonHB«ro  lovr  stack  ooclsts  is 


2.3 


n. 


2.4 


Case  B:   Tp   >  _§_  >  | 


The  stack  height  at  time  T^  is 


SCT^^.  T,  k{^'o,)   . 


2.5 


In  this  case  the  stack  begins  to  decrease  at  a  rate  a  j(  ;  thus  the  timt 
for  the  stack  to  dear  is 

,  2.6 
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The  time  that  a  non-saro  rtaok  ooclsts  la 


Arsrage  low  staek  S. 

G&39  A:         -i-    ^    H.     >    I 


The  average  stack  height  by  definition  is 

sTF)  -.  J~     (^'  set)  dt 


2,7 


.    I 


N  O 


r. 


^    i   °  [T,  Ue-o<)-t<^i2]cit 


T 


P  [T,)l(t-=<^J-tJ?]dt   j  , 


and  after  integratioi  and  slnqplifloatlon 


Z    [g    I-  H    (l-«Ki  "] 


2.6 


Case  B:         -r 

The  arerage  stack  S  is  easily  seen  from  Fig.  2.2  to  be 


2.3  AT»rag«  delay  nT  • 

The  deriratlon  of  the  average  delay  by  the  usual  process  for  com- 
puting averages  Is  rather  umrleldy.  It  Is  nsach  more  conrenlenb  to  com- 
pute T  from  the  ratio  of  the  average  stack  height  S(T)  to  the  average 
landing  rate  J(  (t)  or  the  average  departure  rate  r(t),  all  averages 
being  taken  orer  the  period  T^  during  ifhieh  the  stack  is  nsver  aero. 
That  T  can  very  generally  be  conqpubed  correctly  in  this  manner  may  be 
demonstrated  in  the  follovdng  way. 


SCO 


-t        ti-dt  t, 

stack  Height  vs.  Time  Plot  for  General  Case 

Fig.  2.3 

Consider  a  lovr  stack  ccmosnclng  at  t  ■  0  and  ending  at  t  •  T^     as 
shovm  in  Fig,  2.3.     Assxune  the  variation  of  S(t)  as  arbitrary  bub  subject 
to  the  following  restraints: 

SCO)  -.     S  (T,^=  o, 


S  (t  )    >   O      fofx        o   <  t    <  T", 
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In  terras  of  the  notation  previo««ly  defined,   it  follows  from  Sqs.  2.10 
and  2  Jl  that 

SU)   -.     ]        \n,(.^)  -   ft(t)l      c\± 

^  2J.2 

o  i  -t     6    "T^    . 

This  of  eourse  pots  certain  general  restrictions  upon  r(t)  and  j{  (t), 
since  they  must  lead  to  a  stack  subject  to  Eqs.  2.10  and  2.11.  Granting 
these  general  restrictions,  Eq.  2J.2  is  ralid  regardless  of  the  detailed 
form  of  either  r(t)  or  Jl  (t).  Now  since  S(T  )  -  0, 

rnt)dt  -.     I   a(t)dt  .         ^_^^ 

That  is,  for  such  a  stack  variation,  the  total  number  of  planes  depart- 
ing the  high  stack  equals  the  total  number  landed* 

Consider  the  planes  r(t)dt  which  arrive  in  the  tins  interval  be- 
tween t  and  t  -•■  dt«  They  will  have  to  wait  say,  until  t  ■  t]^  to  land, 
because  of  the  tine  required  to  dear  the  stack  S(t)«  It  follows  there- 
fore that  t]^  is  a  function  of  t  established  by 

which  simply  states  that  all  the  planes  in  the  stack  at  time  t  must 
have  landed  at  a  rate  )(  (t)  by  the  time  t- «  From  Eqs.  2.12  and  2.1if, 
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i  2,15 

0 


Ti^ch  defines  t^  as  a  single  valued  function  of  t.    Thus  as  nd^t  be 
expected,  the  delay  of  a  plane  is  a  function  of  the  tiioe  it  arrives* 
Now  by  definition,  the  average  delay  T  i«  the  total  aggregate 
delay  D  divided  by  the  number  of  planes  delayed,  nhere 

D  '   i^'    (t,-t)  n.it)d±  .- 

0 

That  is,  the  pLsnes  arriving  in  the  time  interval  dt  are  delayed  a  time 
(t^  -  t)  and  the  delays  must  be  sunnned  over  the  period  T    during  vrtiich 
the  stack  is  never  ssero.     Ajod  of  course  the  number  of  planes  delayed  is 
simply  the  total  number  arriving  d\ndng  the  period  T^,     ThuB 


j  V  (tW-t  2a6 


or  in  view  of  Sq,  2.13» 

Now  the  numerator  can  be  integrated  by  parts  as  follows: 
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t)dt  X 


Tc 


trt> 


i  ( 


^(t|dt  dt 


I     • 


2J.8 


By  Eqs.  2.10  and  2,15, 

ti     -    0      vih«n    t  •  0, 

t^     ■    T^    vdMn    t  -  T  . 


2^9 


2.20 


This  slnqply  states  that  the  aircraft  vfhich  arriirs  whsn  the  stack  is  eero 
are  not  delayed  at  all*     Using  Eq.  2J.5  to  eliminate 

o 
ftom  the  last  term  of  Sq*  2*15,  we  hare 


J  'taC*)=lt  =  T^j  ^^t)dt  -J  )f'iit)dijd-t, 


and  by  a  similar  process 


O  _  -         \     o  J 


2.21 


2.22 


Thus 


j(Vt)aWdt=j[|r.CMdtjcl±.J  [j)j(t)dtj 


dt  . 


2«23 


Now  dt]^  is  merely  ftmctionlng  as  a  Tariahle  of  integration  in  the  last 


term  of  Eq.  2.23,  ■<>  ^^at 


—r' 

fU.-t),tt^dt.  Jp;'[.(tw^(t)]dt]d^ . 
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2.24 


Usijig  the  rather  obvloxis  definitions 

s       *  2.25 


*         "  2.26 


S        o 

2.27 


and  Eq8«  2*12  and  2J.3»  v>e  hare  the  desired  alternative  form  of  Eos* 
2J.6  and  2.17: 

■n 


; 


and 


n^ 


)^  U)  n.(t)  2.2s 


Case  A:  J^   i   T^     >  | 

Using  Eqs.  2,4,   2 .26,  and  obserring  Fig.  2  J.,  we  find  that  the  ar- 
erage  landing  rate  is 


iW  :  rr    J"  V«  dt  +   f    ^dt 


s       "^  t;  2.29 
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Thua  from  Eoa,  2.23,  2.29,  and  2.8,  the  average  delay  la 


2e 


2.30 


Case  B:   '  o   ^  ^  > 


From  Flg«  2.2  azid  Eqe.  2.26,  2.2d,  and  2.9»  the  arsrage  delay  la 


n-  - 


.  T 


:^(-^-') 


2.31 


2«4  MfixInniTn  delay  T^   . 

In  order  to  find  the  majcLmum  delay,  the  delay  of  plane  a  arriring 
at  all  possible  times  vrithin  Tg  most  first  be  considered. 


Stack 

Height 

S(t) 


Fig.  2.4  Stack  Height  ts.  Time 
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A8  illustrated  in  Fig.  2.4,  a  group  of  planes  arriving  in  the  time 
interval  (t,  t  +  dt),  with  t  <  T^,  finds  a  stack  of  j!  (€  -  a)t  jUajMa 
waiting  to  land.  Now  the  landing  rate  is  a  X  at  any  time  t  <  Tq,  so 
that  the  delay  for  the  group  of  planes  arriving  during  tims  dt  is 

<   %  '  «<  2.32 

Eq.  2.32  is  only  valid  if  the  delay  'T  <  T^  -  t,  for  if  T  >  T^  -  t,  all 
of  the  planes  will  not  hars  landed  by  T^  and  we  hare  to  consider  the 
change  in  landing  rate  at  T^.  Thus  another  lindting  condition  on  t  in 
Eq«  2«32,  in  additicn  to  the  ccndition 

is  that 


0  ^  1:  i  ^  T,  , 

2-33 


The  real  limiting  ccndition  on  t  in  Eq«  2*32  therefore  depends  upon 
whether 

T^     <   ^   To  2.34 


or 


^  »v.  -     "^     '  o  .  2.35 


3d 


For  the  condition 


T;  >  t  >  i.  X,    (c.se  A) 

^  2.36 

the  planea  arriring  in  the  intenral  dt  will  not  hare  landed  by  T^. 

So  the  number  of  planes  still  due  to  land  before  them  at  T  1« 

Because  these  planes  land  at  a  rate  j(  , 

X 

Now  it  should  be  clear  that  the  maximum  delay  -vdll  occur  for  some 
group  of  pLanes  arriving  at  or  Just  before  Tj.,  where  the  stack  takes 
on  its  maxlBnam  ralue  and  the  landing  rate  remains  small.  The  difficul- 
ties of  finding  the  maximum  near  the  end  point  of  a  linear  e^^uation 
of  the  form  of  'Sr,   2.3^  need  not  be  elaborated  on.  The  mthod  used 
here  is  empirical  and  the  ralue  of  t  obtained  to  maxlmiEe  Er,   2.3^  i^ 

Thus 
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For  the  oondltloo 

the  delay  for  any  plane  in  the  group  is  glren  almply  by  Eo,  2.32,   and 
for  T  to  be  a  maximum  set 

t-V 

Thua 

2.U 
far      Tft      i   ^    >    I  J   Case  B, 

2.5     Additional  delay  of  aircraft  carrier  T„^, 

It  should  be  clear  that  the  extra  tine  needed  by  the  aircraft  car- 
rier to  land  aircraft  v^en  a  low  stack  exists  is  simply 

T. 
where 


ac  ■  '•  "  '^r'  ^'^ 


T      -  N/r. 
r 

Thus  T„^  is  the  additional  time  required  for  the  carrier  to  steam  into 
ao 

the  vdnd  at  high  speeds  axid  should  be  made  a  nrlninwrn  for  the  most  ef- 
ficient utiliBation  of  the  landing  platform. 

Case  A:   -L  >  Zs.  >  I 

From  Ti>is.  2.4  and  2.42, 

T^  (!-•<.)    -  T-n    tl  -  €  )    .  2.43 


ac 
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Case  B:   "^o   >  -L.  >  I 
T"„   ""  o< 

From  EqSo  2o7  and  2.42, 

2c6  Summary  of  results. 

In  the  same  manner,  for  T  <  T  the  calculations  split  up  into 
two  portions,  depending  upon  the  time  at  which  the  stack  finally  clears o 
See  Figs^  2o5  and  2 060  The  complete  results  are  sunmarized  in  Table  2, 
which  gives  the  length  of  the  stack  period,  the  maximum  stack  size,  the 
average  delay,  the  maxiimim  delay,  and  the  aircraft  carrier  delay « 

From  Table  2  the  aircraft  carrier  delay  T^q   is  seen  to  be  in  general 
a  function  of  the  variables  T  ,  T  ,  a  ,  and  £0  It  is  convenient  to 
use  a  normalizing  factor  in  preparing  graphs  for  visual  representation© 
T  was  chosen  as  the  normalizing  factor  for  convenience  because  in  prac- 
tice it  is  usually  fixed 0  A  family  of  curves,  for  a  fixed  value  of  o  , 
showing  T^^^Tj.  vs  the  variable  y  =  Tq/T^.  (to  specify  Tq)  for 
several  values  of  €  as  a  parameter,  were  preparedo  See  FigSo  2o7,  2o8, 
and  2q9o   Other  curves  could  be  prepared  in  a  similar  manner » 

The  general  implications  of  the  numerical  results  may  best  be  ap- 
preciated in  tbrms  of  examples 0  Consider  a  carrier  conducting  landing 
operations  with  the  traffic  density  £  approximately  1,  Now  let  us 
assume  that  during  the  landing  operations  the  landing  rate  has  been 
suddenly  reduced  from  the  normal  rate  in  such  a  manner  that  when  a  =   Oo7, 
Tq  -  0.4  Tj.j  a   -  Oo5,  T^  -  0o2  T,,;  and  a  -  Oo3,  T^  -  0,2  T^.©   In 
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other  words  during  landing  operations  the  landing  rate  has  been  reduced 
in  several  instances  so  that  on  the  average,  say  0o8  of  the  time,  land- 
ings were  being  made  in  an  inefficient  manner  (poor  pilot  technique, 
ine^qperienced  plane  handlers,  etCo)o  Now  from  FigSo  2<,7,  2.8,  and  2o9 
we  find  that  the  additional  delay  to  the  carrier  T^^  to  be  approximately 
0,i*Tj.,  That  is,  the  carrier  must  spend  almost  half  the  time  normally 
taken  to  land  aircraft  steaming  at  high  speed  using  extra  fuel  in  a 
most  inefficient  manners  Taken  over  a  period  of  time  this  could  amount 
to  a  sizeable  wastage,  to  say  nothing  of  the  danger  that  might  be 
present  in  enemy  waters. 

Similarly,  other  values  of  the  parameters  may  be  used  to  compute 
new  delays,  and  in  addition  the  behaviour  of  the  stack  may  be  inves- 
tigated, ioeo,  maximixm  stack,  average  delays,  etc.  These  may  be  impor- 
tant in  determining  fuel  reserves  of  the  aircraft o 

The  hazards  in  vrorking  with  the  traffic  density  £  approximately 
1  have  been  stressed  throughout  the  paper*,  However,  carrier  landing 
operations  are  a  unique  application  of  queueing  theory  and  it  should  be 
clear  that  the  most  efficient  utilization  of  the  landing  platform  will 
only  be  accomplished  with  £  approximately  lo  And  this  can  be  accOTiplifhr. 
ed  only  when  the  pilots  and  other  personnel  are  exercising  the  most  pre- 
cise skill.  Judgment,  etc,  during  the  landing  operationsc  This  is  of 
course  merely  a  statement  of  what  common  sense  would  dictate, 
2o7  Conclusions, 

For  a  simplified  situation  occuring  in  the  landing  of  carrier  air- 
craft, sm  investigation  has  been  made  into  the  variation  of  the  average 
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aircraft  waiting  time,  the  time  wasted  by  the  carrier  during  landing 
operations,  and  the  distribution  of  the  number  of  aircraft  in  the  stack 
for  different  initial  conditions.  The  simplifications  and  assumptions 
necessary  to  arrive  at  a  simple  solution  make  the  problem  somewhat  arti- 
ficial. However,  it  is  not  completely  unreasonable,  because  in  actual 
practice  the  landing  process  may  actually  behave  in  some  instances  in 
a  manner  similar  to  the  process  assumed. 

The  practical  verification  and  application  of  the  results  des- 
cribed in  this  paper  depend  upon  the  specification  of  the  landing  time 
distributions,  the  traffic  density  ^  ,  and  the  minimum  safe  landing 
interval. 

In  comparing  the  present  results  with  actual  carrier  operations, 
many  observations  of  this  process  must  be  made  during  actual  landings 
so  that  the  validity  of  the  assumptions  made  may  be  verified.  In  prac- 
tice, then,  a  great  deal  of  information  is  required  to  check  the 
numerical  model  which  has  been  used  as  a  basis  for  the  result  pre- 
sented here.  If  these  assumptions  are  not  valid  as  based  on  observation- 
al data,  the  model  discussed  in  the  next  secticm  may  be  applicable. 


A7 


2.S  Statiatical  ImpLlcatlona . 

The  difficulties  encountered  in  any  analytical  solutioB  of  a  queue- 
ing  process  with  regular  arrivals  and  a  respectable  distribution  of  ser- 
vice times  is  apparent  from  Lindley»s  analysis  (See  1.6).  In  particuLar, 
if  the  process  is  not  in  statistical  equilibrium,  the  Individual  waiting 
time  distributions  are  difficult  to  obtain  (Eq.  1.19). 

An  alternative  process  is  a  qxiasi-eoperical  Investigation  using  an 
electronic  calculator  as  in  the  LdT  report,  or  if  that  nsthod  is  not 
available,  the  use  of  random  numbers  may  be  substituted.  This  is  a  well 
known  standard  alternative  to  a  purely  theoretical  analysis,  and  its 
basic  principle  is  as  fdlowa.  If  we  read  off  successive  groups  of,  say, 
four  or  five  numbers  from  tables  of  random  numbers,  then  these  can  be 
regarded  as  numbers,  with  four  or  five  decimal  places,  drawn  at  random 
from  a  rectangular  fre-^uency  distribution  with  limits  0  and  1.  Now  all 
continuous  distribution  functions  are  themselves  rectangularly  distribu- 
ted between  0  and  1,  so  that  the  abscissa  of  any  given  distribution 
function  corresponding  to  a  value  of  the  function  taken  fron  a  table  of 
random  numbers  can  be  regarded  as  a  ruantity  drawn  at  random,  with  the 
appropriate  fre'-!uency,  fjpom  the  frequency  distribution  itself.  Thus  for 
a  giiren  distribution  of  landing  time,  for  example  a  gamma  distribution 
with  given  parameters,  n  and  X  say,  we  can  use  random  numbers  to  run 
off  a  series  of  independently  chosen  landing  times.  This  series  can  then 
be  taken  as  the  acttxal  landing  times  of  sxiccessive  aircraft  for  a  hypothe- 
tical aircraft  carrier.  Data  for  a  large  number  of  hypothetical  landing 
patterns  can  be  obtained  In  this  way,  and  the  enqpirical  distributions, 


maana,  «t«indard  deviations,  etc.,  for  rarloua  altenvatiye  landing  aya- 
tema  can  be  worked  out. 

ITlth  this  nethod  of  procedure  it  la  dear  that  we  have  to  decide 
on  the  basic  frenuency  diatribution  of  landing  tinea  beforehand,  though 
T»e  could  consider  a  limited  range  of  type  a  of  diatribution  (rectangular, 
parabolic,  etc.). 

In  order  to  obtain  a  baaic  frequency  diatribution,  obaenrationa  muat 
be  made  of  actuyal  landings  with  the  landing  tixoea  recorded.  The  usual 
curve  fitting  techninuea  can  then  be  employed.  Once  the  proper  diatri- 
butioQ  haa  been  decided  upon  (thia  may  not  be  easy)  the  method  of  maxi- 
BKun  likelihood  may  be  used  to  find  the  unknown  parameters  from  the 
sample  data  (Uood  10  ) . 

For  computational  eaae,  grouping  of  the  frequency  distribution  could 
be  ubilieed.  From  a  practical  point  of  view,  a  truncated  distribution 
may  be  more  realistic  since  the  probability  of  extreme  values  are  ef- 
fectly  eero.  After  the  preliminary  planning,  the  corresponding  dis- 
tribution function  for  each  tjrpe  of  curve  decided  upon  must  be  calcula- 
ted to  four  decimal  places.  Four  figure  numbers  are  then  read  from  the 
table  of  random  numbers,  and  the  value  of  x  (landing  time)  corresponding 
to  the  raid-points  of  the  luterrals  into  which  the  numbers  fell  are  then 
recorded..  In  this  way  aeries  of  landing  times  can  be  obtained  readily, 
simply  by  taking  successive  groups  of  foiir  figure  numbers  from  the 
tables  of  random  numbers  and  writing  down  the  correapondlng  value  of 
X  read  from  the  calculated  table  of  the  appropriate  diatribution  ftinction. 

Now  for  any  aeries  of  landing  tlmea  the  waiting  time  of  each  air- 
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craft  and  the  carrier 'a  idle  or  wasted  time  at  any  time  may  be  calcu- 
lated by  the  usual  set  of  recursive  equations  (see  Fig,  1.1): 

If  ^^,  i«  negative,  that  is  simply  the  carrier's  idle  time.  ..1th  a 
constant  arrival  interval,  the  coirputations  are  simple  but  tedious. 

Finding  the  distribution  of  the  actual  stack  is  a  little  more  com- 
plicated. The  easiest  method  of  computation  seems  to  be  to  consider 
for  how  long  in  the  intervals  (assumed  constant)  between  successive 
arrivals  there  are  0,  1,  2 n  aircraft  waiting.  For  example,  con- 
sider a  constant  departure  interval  of  1  unit.  Consider  five  successive 
aircraft  who  hare  waits  of  0,  0.5,  1,5,  2,  and  1  units.  Assume  there 
is  no  stack  when  the  first  aircraft  enters  the  landing  pattern.  Thus 
the  first  aircraft  doesn't  wait,  so  there  are  no  aircraft  in  the  stack 
and  it  is  empty  for  1  unit,  say.  The  second  aircraft  waits  0,5  units 
of  the  second  interval,  thus  the  stack  is  enqpty  for  the  remaining  0.5 
tmits.  The  third  aircraft  waits  all  of  the  third  Interval  and  0.5  units 
of  the  fourth  interval;  vrtille  the  fourth  aircraft  waits  all  of  the  fourth 
and  fifth  interval;  and  the  fifth  aircraft  waits  all  of  the  fifth  inter- 
val. Thus  it  foUowB  that  In  the  second  interval  there  is  one  aircraft 
in  the  stack  for  1  unit.  Th  the  fourth  interval  there  are  two  then  one 
aircraft  waiting  fer  0.5  units.  It  is  not  too  difficult  to  carry  out 
the  arithmetical  operations  by  using  tables  for  visual  representation 
as  foUows: 
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Waiting  time         Array        Number  of  a/c  In  atack 

y     0    12       2 

1.5  y^y^,^  1   0.5 

2.0  ^\y\  0.5  1 

1.0  y\ 

The  waiting  times  are  grouped  according  to  the  arrival  interral.  The 
times  during  vihich  0,  1,  ,...,  n  aircraft  are  waiting  during  each  inter- 
ral are  shown  in  the  right  hand  table,  and  are  obtained  from  siiccessive 
diagonals  of  the  center  array.  The  figure  appearing  at  the  N.  E.  point 
of  any  diagonal  shows  how  long  the  stack  is  occupiad  by  the  greatest 
number  of  aircraft  in  the  corresponding  interval.  The  latter  nundser  is 
glren  by  the  number  of  nuantities  in  the  diagonal.  The  difference  be- 
tween this  figure  and  the  one  lamediately  to  its  S*  W.  gires  the  length 
of  time  that  one  fewer  aircraft  are  waiting,  etc.  Subtracting  the  figure 
appearing  at  the  S.  W,  point  of  any  diagonal  from  1,0  gives  the  length 
of  time  that  the  stack  is  unoccupied. 

The  above  discussion  erplalns  how  to  obtain  series  of  aircraft  wait- 
ing times,  carrier  idle  time,  and  the  stack  sise.  For  any  particular 
hypothetical  carrier  landing  process  one  can  calculate  the  mean  aircraft 
waiting  time;  the  total  carrier  idle  time;  and  the  time  distribution  of 
the  stack,  ^y  considering  a  series  of  such  proeesses  it  is  then  possible 
to  estimate  what  can  be  expected  to  happen  in  the  long  run,  provided  the 
observational  landing  data  are  valid. 


BIBLIOCKAPHI 


1. 


2.     Kendall,  D.  G, 


3.     Lljudley,   D.  V., 


4*     Adler,   R.   B«  and 
Frlckar,   S.   J., 

5.     B^Ll,   G.  E.   (1949), 


/    6.     Bowen,    S,  G«  and 

F«aro«y,   T.,    (194^) 

7.     Paarcey,  T,,    (194^) 


8.  Feller,   I?,, 

9.  Doob,   J.  L., 
10.     Ikxxl,   A.  M., 


11,  Paley,   R.E.A.C.  and 
V/iener,   N., 

12.  Marshall,    B.  0.,    Jr. 


Operational  Research  '"uarterly  Vol.  3, 
No.  1,   I%pch  1952. 

"SofflB  ptroblems  in  the  theory  of  ^ueues." 
J.  Royal  Statistical  Society,   Series  B, 
Vd.  nil,   No.  2,  1951,   p.  151. 

"The  theory  of  ^ueues  with  a  single  server." 
Proceedings  of  the  Cambridge  PhUosojAiical 
Society,  Vol.  48,    Part  2,   p.   277. 

"The  flow  of  scheduled  air  traffic." 
lOT^FlLE  Reports,   T.  R.,  196,  199. 

"Operational  Research  in  to  air  traffic 
control",    J.  Royal  Aero.  See.  Oct., 
1949,  965. 

"Delays  in  air  traffic  flow",   J.  Royal 
Aero.  Soc.  April,  1948. 

"Delays  in  landing  of  air  traffic", 
J.  Royal  Aero.  Soc.,  Dec,  1948. 

"An  introduction  to  probability  theory 
and  its.  applications".     (1950) 

"Stochastic  Processes".     (1953) 

"Introduction  to  the  theory  of  statistics", 
(1950). 

"Fourier  transforms  in  the  conqoLex 
domain  ."  (1934) 

"OuBusing  theory".  Informal  seminar  in 
Operations  P-esearch  (1952-53),  Seminar 
paper  No.  15. 


52 


R5755 


___L 


J 


R5753 


2S267 

Rippa  

An  applicatioJi  of  queued 
xng  theory  to  carrier  air* 
craft  landing  delays. 


N0236t 
F£  25  63 

FE  2: 


BJilOERY 

^^39 

di>Kay 

^'^'^    70  1 

xQ           7  0  1; 

^^s^^^^zai^ 

=.^1^2-3^9 

12309 

1^^,^ 

25267 
Rip  pa 

An  application  of  queue ing 
theory  to  carrier  aircraft  land- 
ing delays. 


thesR5755missing 

An  application  of  queueing  theory  to  car 


3  2768  000  99077  4 

DUDLEY  KNOX  LIBRARY 


